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A consistent theory to describe the correlated dynamics of quantum mechanical itinerant spins and 
semiclassical local magnetization is given. We consider the itinerant spins as quantum mechanical 
operators, whereas local moments are considered within classical Lagrangian formalism. By appro- 
priately treating fluctuation space spanned by basis functions, including a zero-mode wave function, 
we construct coupled equations of motion for the collective coordinate of the center-of-mass motion 
and the localized zero-mode coordinate perpendicular to the domain wall plane. By solving them, 
we demonstrate that the correlated dynamics is understood through a hierarchy of two time scales: 
Boltzmann relaxation time r e i, when a non-adiabatic part of the spin-transfer torque appears, and 
Gilbert damping time tdw, when adiabatic part comes up. 



Spin torque transfer (STT) process is expected to rev- 
olutionize the performance of memory device due to non- 
volatility and low-power consumption. To promote this 
technology, it is essential to make clear the nature of the 
current- driven domain wall (DW) motionp], [2j]. Recent 
theoretical 0, [5|, @, 0, H, S E| an d experimental [TT| 
studies have disclosed that the STT consists of two vec- 
tors perpendicular to the local magnetization m(x) and 
can be written in general as N = cid x m-\- c^ra x 
d x m\M' The c\ and cvterms respectively come from 
adiabatic [J, 3 and non-adiabatic [5] processes between 
conduction electrons and local magnetization, and the 
terminal velocity of a DW is controlled by not c\ but 
small C2 term. The origin of the C2 term is ascribed to 
the spatial mistracking of spins between conduction elec- 
trons and local magnetization [5]. Behind appearance of 
the C2 term is the so called transverse spin accumula- 
tion (TSA) of itinerant electrons generated by the elec- 
tric current [0, 13] • Now, any consistent theory should ex- 
plain how the adiabatic and non-adiabatic STT come up 
starting with microscopic model. In particular, it should 
be made clear how the TSA caused by the non-adiabatic 
STT eventually leads to translational motion of the whole 
DW. In this letter, to solve this highly debatable prob- 
lem, we propose a consistent theory to describe the cor- 
related dynamics of quantum mechanical itinerant spins 
and semiclassical local magnetization. 

We consider a single head-to-head Neel DW through 
a magnetic nanowire with an easy x axis and a hard 
z axis. Fee electrons travel along the DW axis (x- 
axis). We describe a local spin by a semiclassical vec- 
tor S= Sn=S '(sin cos (/?, sin sin (/?, cos 0) where S = |<S| 
and the polar coordinates and ip are assumed to be 
slowly varying functions of one-dimensional coordinate x 
[Fig. 1(a)]. The DW formation is described by the Hamil- 
tonian (energy per unit area) in the continuum limit, 
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where a is the cubic lattice constant, J is the ferromag- 
netic exchange strength, A = \J JjK and n = y J/K± 



respectively represent the single-ion easy and hard axis 
anisotropics measured in the length dimension. The 
stationary Neel wall (#o = tt/2) is described by no = 
(cos(/?o, sin^o, 0) with ipo(z) = 2 arctan(e x / A ). In the in- 
finite continuum system, the DW configuration has con- 
tinuous degeneracy labeled by the center of mass posi- 
tion, X, of the DW. This degeneracy apparently leads to 
rigid translation of the DW, i.e., no(x) — > no(x — X)[12]. 
As explicitly shown below, however, the translation in 
off-equilibrium accompanies internal deformation of the 
DW. 

The creation operator of a conduction electron is writ- 
ten in a spinor form as c^(x) = (c|(x), c|(x)). By per- 
forming the local gauge transformation c(x) = U(x)c(x) 
with the unitary operator U(x) = e lcr z^o(x)/2 ^ ^ g a 
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FIG. 1: (a) Stationary configuration of local spins (no) asso- 
ciated with a single Neel wall. Labotatory frame x, y, z and 
local frame x, y, z are indicated, (b) Schematic view of the 
transverse spin accumulation (TSA) of itinerant spin s and 
the out-of-plane (0) zero-mode accumulation (OPZA) of local 
spin n. These magnetic accumulations respectively cause the 
non-adiabatic torque T2 and adiabatic torque T±. 
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Pauli matrix) the quantization axis becomes parallel to 
the local spin located at x. Assuming \ad x (po(x)\ ~ 
a/X <C 1, i.e. wall thickness is much larger than atomic 
lattice constant, this procedure leads to the single- 
particle Hamiltonian, 
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where the effective mass of the conduction electron is 
m*. The SU(2) gauge field [1, O is introduced as A z = 

i~ 1 U~ 1 d x U= —(d x (po)a z /2. The conduction electrons 
are assumed to interact with the local spins by a s-d 
coupling represented in the form, 
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Jsd 

a 3 



/CO 
dx s(x) • S(x — X), 
-co 



(3) 



where s and S = Sn are respectively the spins of itin- 
erant and localized electrons. We treat s(x) = ^c^&c as 
fully quantum mechanical operator, while n is a semi- 
classical vector. 

Boltzmann relaxation: let switch on the electric field 
E at t = 0. We introduce the Boltzmann relaxation 
time r e i and the number density of the conduction elec- 
trons fka in the state We assume that the devia- 
tion from equilibrium Fermi-Dirac distribution fo (sua) = 
[exp [(e ka ~ AO /k B T] + l]" 1 is small, where e ka (a =|, j) 
is the single-particle energy, /i is the chemical potential. 
Using standard Boltzmann kinetic equation with relax- 
ation time approximation [6], the distribution function is 
written as 



fka - fo (£ka) + eEr e iv k , 



dfo(Ska) 



(4) 



where the electron charge is — e and the spin-dependent 
velocity is v ka = h~ 1 d£ka/dk. The spin-dependence of 
Ska originates from the SU(2) gauge fields (Az)n anc ^ 
(Az)u- In the process of approaching to stationary cur- 
rent flowing state around the time t ~ r e i, as we will show 
explicitly, the statistical average of the conduction elec- 
tron's spin component perpendicular to the local quan- 
tization axis accumulates and acquires finite value. As 
schematically depicted in Fig. 1(b), this process is ex- 
actly the TSA. The TSA causes an additional magnetic 
field acting on the local spins and exert the non-adiabatic 
torque on the local spins. 

Local spin dynamics: next we formulate dynamics of 
the local spins coupled with the conduction electrons. 
We introduce the SO (x, t) (out-of-plane) and Sep (x, t) 
(in-plane) fluctuations of the local spins around the sta- 
tionary DW configuration rto(x). We say "out-of-plane" 
and "in-plane" with respect to the DW plane. The 
fluctuations are spanned by the orthogonal basis func- 
tions v q and u q as tp (x) = (po (x — X) + Sep (x — X) and 
(x) = it/2 + SO (x - X) , where 

pOO POO 

Sip (x) = / dq r]q(t)Vq (x) , SO (x) = / dq £ q (t)u q (x) . 



(5) 

At this stage, X is not a dynamical variable, but just 
a parameter. The basis functions obey the Schrodinger 



equations, (JS 2 /2) (-<9 2 -2X~ 2 sin 2 (^o +A~ 2 ) v q (x) = 
e*u q {x) and (JS 2 /2) (-d 2 -2A~ 2 sin 2 y?o +A~ 2 +k~ 2 ) 
u q (x) = e e q u q (x) . Both and cp modes consist of a sin- 
gle bound state (zero mode) and continuum states (spin- 
wave modes). The dimensionless zero mode wave func- 
tions are given by uq (x) = vq (x) = &o(x), where 



(6) 



®o(x) =\l —d x (pv(x) = 



2A cosh(x/A), 



with the corresponding energies respectively given by 
£ e = JS 2 /(2k 2 ) and Sq = 0. The normalization is 

given by a -1 j^dx [$>o(x)] 2 = 1. Although to excite 
the out-of-plane (0) zero mode costs finite energy gap 
a® coming from the hard-axis anisotropy, we still call 
this "zero mode." The spin- wave states have energy 
dispersions given by e e = \ JS 2 (q 2 + A -2 + k~ 2 ) and 

£ q = ^JS 2 {q 2 + A -2 )- Because the zero mode and the 
spin-wave states are orthogonal to each other and sepa- 
rated by the anisotropy gaps, the spin- wave modes are 
totally irrelevant to a low energy effective theory. There- 
fore, we ignore the spin-wave modes from now on. 

Out-of-plane zero-mode(OPZ) coordinate ^o- in order 
to obtain the correct form of the dynamical Hamiltonian, 
one has to regard the variable X as a dynamical variable 
X(t) and replace the zero mode coordinate r]o with X(t). 
Following this idea, the zero-mode fluctuations should be 
given by, 



<p(x,t) = (po [x-X(t)\ , 
0(x,t) =7r/2 + fo(*)*o [x- 



X(t)}. 
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Eq. ((HJ) is a key ingredient of this letter, which has 
never been explicitly treated so far[T4[. That is to say, we 
naturally include the out-of-plane(OPZ) zero-mode, in 
addition to the in-plane (cp) zero-mode replaced by X (t). 
The zero- mode wave function $0 [x — X (t)] serves as the 
basis function of the ^-fluctuations localized around the 
center of the DW and £o(£) is the OPZ coordinate. Now, 
our effective theory is fully described by two dynamical 
variables X(t) and £o(£) which naturally give physical 
coordinates along the Hilbert space of orthogonal and 
cp fluctuations. As we will see, we have £o(£) 7^ only 
for inequilibrium current flowing state under E ^ [Fig. 
2(a)]. 

It is here important to note an essential difference be- 
tween Tatara and Khono's approach^ and ours. Tatara 
and Khono used X(t) and the weighted average, 0o(t) 
— J^^dx 0(x, t) sin 2 cp[x — X(t)], as dynamical vari- 
ables. Later, they systematically used complex coordi- 
nate £ = e^tan(0/2) and described the fluctuations in 
the form £ = e -<x,t)+i^ +v[x-x(t)] [| (their notation is 
reproduced by putting ^ tt/2 — 0, (p — > ip in our nota- 
tion). In our understanding, these descriptions inevitably 
cause redundant coupling between u and v modes in Eq. 
([5]). Actually, our natural choice of the dynamical vari- 
ables is essential to appropriately derive relaxational dy- 
namics described b y the following equations of motion 
given by (fT2a) and (|12b|) . 

Equations of motion of the DW: now, we construct 
an effective Lagrangian C= £dw + £sd to describe the 
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FIG. 2: (a) Spatial profile of the polar angles <p(x,i) = 
(po[x - X (t)} and 0(x,t) = tt/2 + £o(*)$o [x - X (t)} in the 
current flowing state, (b) Linear dependence of £o and X (t) 
on the current density j. (c) Single- particle propagation (rep- 
resented by solid line) with spin flip process by the s-d inter- 
action (represented by wavy line) which leads to the STT. 



DW motion and resultant equations of motion (EOM). 
Using ([7j) and (jSj), the local spin counterpart is given by 
£dw = f^dx (cos 6 — 1) (p—Hr>w explicitly written 

as 
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To understand the effect of the s-d coupling, it is use- 
ful to note n[6o + 50, (po + 5cp] ~ no —e z 56 —no50 2 /2, 
where we dropped Sep because this degree of freedom is 
eliminated by the global gauge fixing [lj- We have thus 
s-d Lagrangian, 

C sd = a- 3 J sd S (To ~ <S||£ 2 /2) , (10) 

where, T§ [X (t)] = j^^dx n [x — X (t)] • (s(x,t)) 
and 5|| [X (t)} = J^dx {$ [x - X (t)}} 2 n [x-X (t)} • 
(s(x,t)). Finally, to take account of dissipative 
dynamics, we use the Rayleigh dissipation function 
W Ray iei g h = § S-oodx n 2 explicitly written as 



a hS . - 

WRayleigh = ~ ~^ I < + J X ' 



(11) 



where a is the Gilbert damping parameter. It is simple 
to write down the Euler-Lagrange- Rayleigh equations, 
d(dL/dqi) jdt — dL/dqi = —dW/dqi, for the dynami- 
cal variables q\ = X and q2 = £o- We obtain the EOMs 
which contain the dynamical variables in linear order, 
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(12a) 
(12b) 

(13a) 
(13b) 



respectively give the non-adiabatic STT and longitudinal 
spin accumulation 0. The statistical average of the con- 
duction electron's spin component is denoted by (•••). 
The gauge-transformed spin variables are introduced by 
s(x) = U^lx-Xit)} s(x)U[x-X(t)\ which has lo- 
cal quantization axis tied to the l ocal spin at the po- 
sition of x — X(t). To obtain Eq. (|13ap . we used rela- 
tions d x \\Q [x — X(t)] = —d x (po(x) e z x no [x — X(t)] and 
(s y ) = — (s x ) smcpo + (s y ) cos(/?o- The relation (|13aj) im- 
plies that the translation of the DW (x — > x — X) natu- 
rally gives rise to the TSA, (s y ), along the local y axis. 
The appearance of (s y ) causes local magnetic moment 
which triggers the local spins to precess around the local 
y axis and consequently produce finite deviation of the 
polar angle 50 = — 0q. It is seen that upon switching the 
external electric field, the deviation 50 relaxes to finite 
magnitude in the stationary current- flowing state, i.e., 
the OPZ coordinate £o(t) accumulates and reaches finite 
terminal value £q . We call this process out-of-plane zero- 
mode accumulation(OPZA)as schematically depicted in 
Fig. 1(b). This effect is physically interpreted as appear- 
ance of demagnetization field phenomenologically intro- 
duced by Doring, Kittel, Becker [15|. and SlonczewskipQ]. 
It is also to be noted that we ignored the term dS\\/dX. 
This simplification is legitimate for the case of of small 
sd-coupling. 

G ilbert relaxation: coupled equations of motion ()12aj) 
and (|12b|) are readily solved to give relaxational solutions, 



Co = £o(l- e " t/rDW )> V = X = V*(l-e- t/TDW ), (14) 
where the OPZA reaches the terminal value, 



_ 1 v/oA/2 J sd T ± 
so — ~~~ 
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and correspondingly the terminal velocity of the DW 

reaches V* = zJsdXL- The relaxation time of the 

DW magnetization, tdw? is given by 



k 2 a?JS + JsdS\\ 
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This result clearly shows that the DW magnetization try 
to relax through the Gilbert damping toward the direc- 
tion of the newly established precession axis. We s tress 
that without the OPZ coordinate £o in Eqs. (|12aj) and 
(|12bp . only the terminal velocity is available and the tran- 
sient relaxational dynamics is totally lost. 

As depicted in Fig.2(a), the OPZA [Eq. JTJJJ)] gives rise 
to finite out-of-plane (z) component of the local spin, 



n z (x, t) = COS0 : 



1 (K\ 2 J sd 



2a 
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1 



tT±. 



J JS cosh[(x - X(t))/X\ 

(17) 

The resultant local spin S± = Se z n z (x,t) gives the de- 
magnetization field phenomenologically treated by Slon- 
czewski and gives rise to the adiabatic torque T\ = 
cid x n(x) = c\ (d x (po) (— sin cpo, cos cpo, 0). At the in- 
terface of the DW boundary, (po = tt/2 and T\ = 
ci (d x (po) (—1,0,0), i.e., the adiabatic torque rotate the 
local spin to counterclockwise direction when the electric 
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current flows in the (1, 0, 0)-direction. As is clear from 
the above discussion, this adiabatic torque is established 
after the stationary current-flowing [j = (ne 2 T e \/m*)E) 
state establishes the non-adiabatic torque, 71. Around 
the time scale of t ~ T e i+TDw? the whole system (in- 
cluding conduction electrons and DW) reaches non- 
equilibrium but stationary state. In this state, the DW 
magnetizations continuously feel the OPZA and macro- 
scopically rotate around it. This process exactly corre- 
sponds to stationary translation of the DW. 

Computation of 71: the final step is to compute an 
explicit form of 71. By taking Fourier transform c~ka(t) 
= -^j^ k e %kx £o{%^\ and retaining only the momentum 
conserving process, we have 

f7i/a 



T±.=\\ dk ReG< M (t,t), (18a) 

(18b) 



t/q 



Si 



dk lmG< M (t,t). 

Z7r J-ir/a 



Here, the expectation values are computed by using the 
lesser component of the path-oriented Green function 
G ^,fcV'(M') = ify'A^kait)), where t (f) is defined 
on the upper (lower) branch of Keldysh contour. Since 
S\\ does not play an essential role, we pay attention to an 
essential quantity 71. To evaluate the Green functions, 
we perturbatively treat the s-d coupling and write down 
the Dyson equation. Then, we truncate the Dyson equa- 
tion by using the Born approximation including the s-d 
coupling in linear order which causes a single spin flip 
process [Fig. 2(c)] and gives rise to off-diagonal compo- 
nent in spin space, 



(M) 



2 £/e T - Ski ~ ^ 
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To obtain the explicit form of Ska, we write the single- 
particle Hamiltonian ([2j) in Fourier space and obtain 
T~Le\ = 7~to +7^gauge> where Ho represents free conduction 
and 7Y g auge comes form the second term in Eq. (j2j). By 
retaining only momentum conserving process, we have 
^ei = J2k,a £ kacla^ka, where s khi = h 2 (k =F 5k) 2 /2m* , 
where the shift of the Fermi wave numbers due to the 
background DW is given by 5k = tt/ (2a). 

Using Eqs. (|4]), (O, and (|18a|) . we finally obtain the 
STT which points in the z-direction, T\ = 71 e z , where 
its magnitude is given in a form, 



T, = 



1 Js, 



1 



4 k B T cosh 2 [(s - /i) /2k B T] jo 
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where jo = ^neh/ (nam*) and eo = h ir /(8m*a ) corre- 
sponds to the chemical potential at half-filling. We have a 
master formula which gives relation between the current 
density and the terminal velocity of the DW, 
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8a h k B T cosh 2 [(s - /i) /2k B T] jo 
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As shown in Fig. 2(b), we see there is no threshold for the 
velocity, which is consistent with the result obtained by 



Thiaville et al.[10]. Standard choice of parameters, jo — 
10 16 [A-m- 2 ], A = 10- 8 [m], a = 10~ 2 , j ~ 10 n [A- m - 2 ] 

give a rough estimate V* ~ -100 (J s d/k B T) 2 [m/s]. Of 
course, to pursuit more quantitative result needs numer- 
ical estimation of 71 taking account of real band struc- 
ture. 

It is essential that the Gilbert damping coefficient, a, 
enters Eq. ([2T]) . The relaxation process of the DW dynam- 
ics is governed by the Boltzmann relaxation followed by 
the Gilbert damping in hierarchical manner. As summa- 
rized in Figs. 2 (a) and (b), in our treatment, it is crucial 
to recognize that the OPZ coordinate £o acquires finite 
value (i.e., accumulation) only for the current flowing 
state which is non-equilibrium but stationary. This is the 
case where dynamical relaxation leads to finite accumula- 
tion of physical quantities which are zero in equilibrium. 
Although essential role of the sliding mode to describe 
localized spin dynamics was pointed out before [1, [l2| 
and importance of out-of-plane canting of the local spins 
was stressed[l, 8], the OPZA presented in this letter has 
not been discussed before. For example, the sliding mo- 
tion in Ref.[l2| does not contain internal deformation 
of the DW. The OPZA is an outcome of time-reversal- 
symmetry breaking by electric current. This interpre- 
tation seems natural because current-flowing state is off 
equilibrium. 
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